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In liquid mixtures and other binary systems at low temperatures the pure phases 
may coexist, separated by an interface. The interface tension vanishes according 
to a = <7o(l — T/T c )^ as the temperature T approaches the critical point from 
below. Similarly the correlation length diverges as £ = /-(l — T/T c )~ u in the 
low temperature region. For three-dimensional systems the dimensionless product 
R- = <7q/£ is universal. We calculate its value in the framework of field theory in 
d = 3 dimensions by means of a saddle-point expansion around the kink solution 
including two- loop corrections. The result i?_ = 0.1065(9), where the error is mainly 
due to the uncertainty in the renormalized coupling constant, is compatible with 
experimental data and Monte Carlo calculations. 
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For a statistical system near a critical point various measurable quantities X obey a 
singular behaviour as a function of the temperature T like 



with a critical exponent e, where 



and T c is the critical temperature. For a given observable X the critical exponent is a 
universal quantity and assumes the same value for systems belonging to the same univer- 
sality class. The critical amplitude X , however, is not universal and varies from system 
to system, depending on the microscopic details of the Hamiltonian. 
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The critical exponents of a universality class are not independent of each other but 
obey a number of scaling and hyperscaling relations. For example, in the low temperature 
region the exponents /3, 7 and v belonging to the magnetization 



the susceptibility 

and the correlation length 

are related by 



M ~ Bt p , (3) 
X ~ C_t~\ (4) 

£ ~ S-t~ v (5) 



2(3 + 7 - dv = 0, (6) 
where d = 3 is the number of dimensions. Therefore in the combination 

the exponents cancel and ur is expected to approach a finite value 

^ (8) 

when the critical point is approached from below. It is a another consequence of the 
scaling hypothesis that such combinations of critical amplitudes are also universal |], ||. 
They are generally called amplitude ratios. For a review see 

In the last decades much interest has focussed on critical indices, whose values are 
known from various methods rather accurately by now. From a phenomenological point 
of view the amplitude ratios are, however, at least as interesting as the indices. They are 
well accessible experimentally and their numerical values are often more characteristic for 
the universality classes as the variation between different classes are larger. 

In this article we consider an amplitude ratio related to the interface tension in the 
universality class of the three-dimensional Ising model. In various binary systems at 
temperatures T below T c interfaces (domain walls) may be present, separating coexisting 
phases. The interface tension r is the free energy per unit area of interfaces. As T 
increases towards T c the reduced interface tension 

° - h> < 9 » 

where k is Boltzmann's constant, vanishes according to the scaling law 

or ~ <7 t M . (10) 

Widom's scaling law ||, |], 

fi = 2u, (11) 
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relates the universal critical exponent \x to the critical exponent v of the correlation length 
£. Associated with this law is the universal dimensionless product of critical amplitudes 

R- = <?ofl. (12) 

In this article we consider the correlation length as defined by means of the second 
moment of the correlation function || , in contrast to the "true" (exponential) correlation 
length. Numerically they differ by less than 2 percent 0. 

The amplitude a has been studied experimentally (see 0, as well as theoretically 
The presently most accurate Monte Carlo result has been obtained by Hasenbusch and 
Pinn §. 

In the theoretical toolbox we also have the three-dimensional Euclidean 4 -theory, 
which is believed to be in the same universality class as binary systems and the Ising 
model. Therefore it should describe the universal properties of these systems correctly. 
The scalar field 4>(x) represents the local order parameter which in the case of binary fluid 
mixtures is proportional to the difference of the concentrations of the two fluids. 

On the classical level an interface in a system of cylindrical geometry is represented 
by a classical solution of the field equations. It is a saddle-point of the Hamiltonian. 
Taking thermal fluctuations into account amounts to performing an expansion around 
the saddle-point in the functional integral. 

In the field theoretical framework the interface tension, in particular the universal 



ratio has been investigated by means of the e-expansion in 4 — e dimensions [10|, 
II| and directly in 3 dimensions |TI|. In both approaches the calculations were done 



in the one-loop approximation (see below). Whereas the results from the e-expansion 
are afflicted by convergence problems and show large deviations from experimental and 
Monte Carlo results, the three-dimensional field theory leads to relatively small one-loop 
corrections and more reasonable numbers. But higher-loop corrections may spoil this 
situation. Therefore, in order to get a better impression of the numerical convergence and 
to obtain more precise estimates it is highly desirable to know the two-loop contribution 
to R-. In this article we present the result of a two-loop calculation of the universal 
amplitude ratio R- in the framework of three-dimensional 4 -theory. 

The Hamiltonian Ti,, which is called action in the context of Euclidean field theory, in 
the broken symmetry phase is written in terms of the bare field 0o as 

H = Jcd 3 x, £ = ^odxfo + V(<p ), (13) 
where the double-well potential 



has its minima at 



±^o = ±a (15) 
\ 9o 



The parameters are defined such that the value of the potential at its minima is zero and 
mo is the bare mass. The renormalized mass 



m R = 1/f 



(16) 



is the inverse of the second moment correlation length. It is defined together with the wave 
function renormalization Zr through the small momentum behaviour of the propagator: 



1 



JR 



{ m 2 R + P 2 + 0(p 4 )}. 



The renormalized vacuum expectation value of the field is 



Vr 



J R 



-1/2, 



(17) 



(18) 



where v is the expectation value of the field 
coupling we adopt the definition 



For the dimensionless renormalized 



ur 



9r 



mR 

i > 

V R 



(19) 



which in the language of statistical mechanics corresponds to Eq. (^). 

The basic idea behind the calculation of the interface tension is its relation to the 



energy splitting due to tunneling in a finite volume. We refer to |13], for details. In 
a rectangular box with cross-section L 2 the degeneracy of the groundstate of the transfer 
matrix is lifted by an "energy" splitting E 0a . This gap depends on L according to 



Ei 



(20) 



= Cexp{-aL 2 } 

where a is the (reduced) interface tension fT4] , |1"5| , IB] , |I~2 

The energy splitting E§ a can be calculated in a semiclassical approximation, which 
amounts to a saddle-point expansion around the classical kink solution 



(f) c (x) 



3m 2 m 

a tanh — [xs — a) 

\ 9o 2 



(21) 



of 4 -theory, where a is a free parameter specifying the location of the kink. The classical 
energy of a kink is 

(22) 



m: 



n c = 2 ^L 2 . 

9o 

The kink interpolates between the two field values at the minima of the potential and 
represents an interface separating regions with different local mean values of the field. 
In the two-loop approximation the functional integral 



-H[<h] 



v<p 



(23) 
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(a factor of 1/kT has been absorbed in Tl) with appropriate boundary conditions is 
calculated by expanding the energy TC[4>] around the kink solution C up to order g and 



evaluating the integral by the saddle point method. For details see fll7| . An analogous 
calculation for the case of the anharmonic oscillator in quantum mechanics has been 
performed in [R]. 

The zero mode belonging to the translations of the center of the kink, i.e. to shifts in 
a is treated by the method of collective coordinates and leads to a nontrivial Jacobian J. 
For the energy splitting one obtains 



E, 



H c /det'M\- 1 / 2 



Oa 



2vr VdetM, 



x exp < -H c + 



+ 




+ 



CO - CO 




1 
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+ 0{gl) , 



where 



M 



-d^ + m„- ^m 2 cosh 2 (^r^3 



is the operator of quadratic fluctuations around C and 

M = -d^ + ml 



(24) 



(25) 



(26) 



The prime in det' indicates a determinant without zero-modes. The two-loop contributions 
are displayed as Feynman diagrams. The propagators are 



= G(x,y), 

= G (x,y), 

where the propagator in the kink background 

G(x,y) = (x\(Mr 1 \y) 
is the Greens function of M without zero mode and 



G (x,y) = (x\M 



-ii 



(27) 



(28) 



(29) 



is the usual scalar propagator. Remember, however, that the propagators refer to a system 
with finite cross-section L 2 . 
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The vertices are 




V - - go, 



® - ~TtMx). (30) 

ric 

The last vertex comes from the Jacobian J. 

The spectrum of the fluctuation operator M is known exactly. Owing to this the 
determinants can be evaluated analytically with the help of heat kernel and zeta-function 
techniques [T2|]. They yield the one- loop contribution to the interface tension. 



Much more involved is the calculation of the two-loop contributions. Although we have 
an expression for G(x, y) (covering one and a half page) it turned out to be advantageous 
to use the Schwinger representation 



/•oo 

G= / e- tM dt (31) 
Jo 



and to write the kernel exp(— tM') in the spectral representation. The calculations have 
been done analytically as far as possible. In the later stages some infinite sums and low- 
dimensional integrations have been done numerically. Details of the calculation can be 



found in [17 



The most difficult piece, of course, was the true two-loop diagram |CJ) [ ■ ^ con tains 
ultraviolet divergencies, which we isolated by means of dimensional regularization in d — 
3— e dimensions. After some tedious calculations (we warn the curious reader) we obtained 
the two-loop contribution as a function of L up to terms of order L°. Whereas individual 
diagrams produce terms of the form L 2 (log(moL)) 2 and L 2 \og(m L), they cancel in the 
total sum. The leading L-dependence is then proportional to L 2 as is required by the 
finiteness of the interface tension. 

The ultraviolet divergence is removed by renormalization as usual. Expressing the 
bare parameters m and #o m terms of their renormalized counterparts m# and in the 
two-loop approximation indeed cancels the pole in 1/e. The final result for the interface 
tension at L = oo is 

-^HsME)'^}- (32) 



with 

1 / 3 \ 37 
° 11 = I ( 3 + I log 3 J ~ 32 = -°- 2002602 ( 33 ) 

and 

<r 2i = -0.0076(8). (34) 
Using (p~6|) the desired amplitude ratio R_ is obtained by evaluating the function 

f(u R )=a/m R (35) 

at the fixed point value u R = u* R , i.e. 



R- = f(u R ) = — l + ^-f + aa -f +0«) . (36) 



Mjj | 47T \47T 



The most recent results for tt|j are 



u* R = 14.3(1) (37) 



from Monte Carlo calculations @ and u* R ~ 14.2 from three-dimensional field theory 
| i~9f . Earlier estimates were 14.73(14) from low-temperature series |20] and 15.1(1.3) 



used in [0 . For these numbers the two- loop contribution to R_ is about 1 % while 



the one-loop contribution is about 24%. Although the apparent numerical convergence 
is surprisingly good we have also applied Pade and Pade-Borel approximations to the 
quadratic polynomial appearing in f(un) in order to get The dependence of R on 
u* R is illustrated in Fig. 1, where the average of the three Pade approximants and the 
corresponding error band is displayed. 

Table 1 shows the results of the Pade and Pade-Borel approximations evaluated at 
three of the values for u* R mentioned above. The quoted errors are due to the error of a%. 









/ 




/[0,2] 


/[1,1],PB 


/[0,2],pb 


15.1 


0.0991(2) 


0.0989(2) 


0.1011(1) 


0.0991(3) 


0.10454(5) 


14.73 


0.1025(2) 


0.1024(2) 


0.1044(1) 


0.1025(3) 


0.10774(5) 


14.3 


0.1066(2) 


0.1064(2) 


0.1084(1) 


0.1066(3) 


0.11166(5) 



Table 1: Results for R_ from Pade and Pade-Borel approximations 
evaluated at three values for u* R . 



The average of all approximants at u* R = 14.3 is R- = 0.108(2). The [0,2] approx- 
imants appear to be off the rest. Leaving them out yields R- = 0.1065(1). Taking the 



error of u* R into account we obtain 



i?_ = 0.1065(9). (3* 
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Figure 1: R- as a function w^. The solid line is the average of the three 
Pade approximants. The thin lines indicate an error estimate of 0.001. 

Since we do not know the size of higher-loop contributions the quoted error mainly reflects 
the spread due to the uncertainty of u* R . 

For comparison the Monte Carlo calculations of Hasenbusch and Pinn || yield 

= 0.1040(8). (39) 

In view of the remarks made above we consider the results as being compatible. 

Experimentally the universal amplitude combination R + = o"o/+ ; where /+ is the 
amplitude of the correlation length in the high temperature phase, has been measured for 
various binary systems, see |§. In order to compare with i?_ the universal ratio /+//_ 
has to be employed. Using /+//_ = 1.95(2) from recent Monte Carlo calculations [[7] or 
/+//_ = 1.99(2) from field theory (see |19| and the remark in the conclusions of 0) we 
obtain for R + the numbers 0.40(1) and 0.42(1), respectively. This compares well with 
the recent experimental result of 0.41(4) for the classical cyclohexane-aniline mixture 0. 
Previous experimental results are summarized in 0] as 0.37(3). 
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